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Abstract 

We show that the dilepton spectra measured by the CERES collaboration in Pb 
■ + Au interactions for various charge multiplicities can be reproduced by a hadronic 

\ initial state with reduction in the masses of the vector mesons in the thermal bath. 

. Though such an effect may also be achieved by a large broadening of the spectral 

function we show that the photon spectra is insensitive to this. It is found that for 
higher multiplicities a good description of the data can also be obtained with quark 
gluon plasma initial state if the reduction of the vector meson masses in the mixed 
and hadronic phases is taken into account. We observe that a thermal source with 
initial temperature ~ 200 MeV can reproduce the observed enhancement in the low 
mass region of the dilepton spectra. It is not possible to state which one of the two 
' initial states (QGP or hadronic) is compatible with the data. These findings are in 

agreement with our earlier results obtained from the analysis of the WA98 photon 
| spectra. We estimate the number of n — n collisions near the p-peak of the dilepton 

. spectra and argue that thermal equilibrium may have been achieved in the system, 

justifying the use of hydrodynamical model to describe the space time evolution. 



o 

1 Introduction 

One of the basic aims of research in highly relativistic nuclear collisions has been the 
creation and detection of Quark Gluon Plasma (QGP). Such an endeavour started at the 
BNL AGS, continued at the CERN SPS and is presently being pursued at the heavy ion 
experiments at the BNL RHIC. The principal motivation behind such a sustained and 
organized effort is that apart from being a confirmation of a long standing prediction 
^ ■ of thermal QCD, the production of deconfined matter would actually recreate in the 

laboratory conditions which existed in the microsecond universe. A variety of signals 
of this exotic state of matter have been proposed and studied extensively over the last 
two decades. The suppression of the J particle, the enhanced production of strange 
particles, specially strange antibaryons, excess production of photons and lepton pairs 
and the formation of disoriented chiral condensates are some of the well accepted probes 
of a deconfinement and/or a chiral symmetry restoring phase transition. In the recently 
concluded Pb run at the CERN SPS some of these, specially the magnitude of J 
suppression has been claimed to be inexplicable in a non-QGP scenario Q. 

Owing to their long mean free paths, electromagnetic signals viz. photons and dilep- 
tons are by far the most direct probes of the collision ||, ||, ||, ||, ||. In fact they possess 
the unique property of being sufficiently less prone to secondary interactions because of 
their electromagnetic nature and yet are strong enough to be experimentally detectable. 
Moreover, unlike hadronic signals, photons and dileptons are produced all through the 
evolution process. In this connection one must bear in mind that even if QGP is produced 
in the initial stages, the quarks and gluons will eventually form colour singlet hadrons 
before traveling to the detectors. Hence, in order to unmask the emissions from the QGP 
one must have a rather accurate estimation of the contributions from hadronic sources. 
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In this aspect, the dilepton spectrum enjoys a clear advantage. This is because, though 
both photons and dileptons couple to hadrons through spin one "1 mesons, the 

dilepton spectra exhibits a resonant structure which, in the low mass regime includes 
the p and the u mesons. Now, it has been emphasized that the properties of the vector 
mesons undergo nontrivial medium modifications in a hot and /or dense medium such as 
likely to be produced in relativistic nuclear collisions. Consequently the spectral modi- 
fications of these vector mesons would be clearly revealed in the invariant mass spectra 
of the dileptons through the shift of the resonant peaks. 

Let us recall the various sources which produce dileptons in relativistic nuclear col- 
lisions. The principal source of thermal dileptons from the QGP and hadronic phases 
are the quark-antiquark and pion annihilation processes respectively. In the high mass 
region these compete with Drell-Yan pairs which are produced in primary interactions 
between incoming partons in the very early stages of the collision. The J/^> peak at 
3 GeV marks the cut-off scale for thermal pairs from the plasma. Around this region 
the decays of D and B mesons also become an important source 0] . The vector mesons 
which decay both during the expansion and after freeze-out can be identified from their 
characteristic peaks in the spectrum. Below the vector mesons, Dalitz decays of ir°, n, 
rjf and uj mesons provide the dominant source for dilepton production. 

The observed enhancement of low-mass dileptons in Pb+Au collisions at the CERN 
SPS compared to the yield from hadronic decays at freeze-out as reported by the CERES 
Collaboration [|| has triggered a host of theoretical activities. Though the pion annihi- 
lation channel 7r + 7r~ — > l + l~ accounts for a large fraction of this enhancement, it turns 
out that a quantitative explanation of the data requires the incorporation of medium 
modifications of the vector mesons. The enhancement of the dilepton yield in the low 
invariant mass region was well reproduced by using reduced vector meson masses in the 
medium [H 10 1. Rapp et al have used a large broadening of the p meson spectral 



function due to scattering off baryons in order to explain the data. However, both ap- 
proaches rely on a high baryon density for the dropping mass or the enlarged width 
of the p meson. Koch et al [|l^] finds very little effect due to baryons on the dilepton 
yield in the kinematic range of CERES experiment. It is worth emphasizing here that 
as yet it has not been possible to explain the observed low-mass enhancement of dilep- 
tons measured in the Pb+Au collisions [|R| as well as in the S-Au [14] collisions at the 
CERN SPS in a scenario which does not incorporate in-medium effects on the vector 
meson mass. In short, the CERES dilepton data indicate a substantial change in the 
in-medium hadronic spectral function, either mass shift or broadening but fails to make 
any distinction between them. 

We have shown earlier [15] that the WA98 photon spectra can be explained by 
either of the two scenarios of relativistic nuclear collision: (a) A +A -^QGP— >Mixed 
Phase— ^Hadronic Phase or (b) A + A -^Hadronic Matter; with downward shift of vector 
meson masses and initial temperature ~ 200 MeV. The photon yield is insensitive to 
the broadening of vector mesons [16], although the CERES dilepton data admits such 
a scenario as mentioned above. In the present work we would like to see whether the 
CERES dilepton data can be described by the same scenarios ( a and b) and with similar 
initial conditions which reproduce the WA98 photon spectra. The effects of the thermal 
shift of the hadronic spectral functions on both photon and dilepton emission have been 
considered in Ref. (l(| for an exhaustive set of models. However, an appreciable change 
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in the space-time integrated yield of electromagnetic probes was observed for universal 
scaling and Quantum Hadrodynamic (QHD) model. Therefore, in the present calculation 
we will consider these two cases. 

The paper is organized as follows. In the following section we will briefly outline 
the phenomenology of medium effects on the vector mesons in the thermal environment 
which we have considered. In section 3 we will consider the static dilepton rates due 
to different processes. Then in section 4 we will describe the dynamics of space-time 
evolution followed by the results of our calculation in section 5. We will conclude with 
a summary and discussions in section 6. 



2 Medium Effects 

A substantial amount of literature has been devoted to the issue of temperature and/or 



density dependence of hadrons within various models (see |16, 0' @ Hp 13 f° r a 
review). In Ref. [|l6| we have discussed the effects of spectral changes of hadrons on the 
electromagnetic probes in detail. In this work we will consider medium modifications of 
vector mesons in two different scenarios where non-trivial effects on the dilepton spectra 
were seen. These are: the universal scaling hypothesis and the Quantum Hadrodynamic 
(QHD) model. 

In the scaling hypothesis, the parametrization of in-medium quantities (denoted by 
*) at finite temperature, T and baryon density, n# pi] is 

^^4_o. 2 -Vi-^\ (1) 

my f V U \ n B J \ T c J 

where V stands for vector mesons, fv is the coupling between the electromagnetic current 
and the vector meson field, ujq is the continuum threshold, T c is the critical temperature 
and n'g is the baryon density of normal nuclear matter. Mass of the nucleon also varies 
with temperature according to Eq. (|l]) (pseudo scalar masses remain unchanged). The 
particular exponent A = 1/2 (1/6) is known as Nambu (Brown-Rho) scaling. We will 
use A = 1/2 in our calculations. Note that there is no definite reason to believe that 
all the in-medium dynamical quantities are dictated by a single exponent A. This is the 
simplest possible ansatz (see |16| for a discussion). The effective mass of a\ is estimated 
by using Weinberg's sum rules 



In the Quantum Hadrodynamic model [23, 24 1 of nuclear matter the vector meson 



properties are modified due to coupling with nucleonic excitations. The nucleons interact 
through the exchange of scalar a and the vector to mesons and their mass is modified due 
to the scalar condensate. This is evaluated in the Relativistic Hartree Approximation 
(RHA). Coupling with these modified nucleonic excitations induce changes in the p 
and ijj meson masses. This modification is contained in the meson self energy which 
appears in the Dyson-Schwinger equation for the effective propagator in the medium. 
The interaction vertices are provided by the Lagrangian 

£%n = 9VNN (n^NV* - ^-Na, v r a Nd^ , (2) 

where VJf = {uj^, p 5 "}, Mn is the free nucleon mass, N is the nucleon field and r a = {1, r}, 
t being the Pauli matrices. The real part of the vector meson self-energy due to NN 
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polarization is responsible for the mass shift. The details of the calculations can be 
found in our previous work 2J, |26j (see also p7j ) and we do not reproduce them here. 



3 Dilepton Emission Rate 

In this section, we briefly recapitulate the main equations relevant for evaluating dilepton 
emissions from a thermal source. The emission rate of dileptons with four momentum 
p^ can be expressed in terms of the imaginary part of the retarded current-current 
correlation function as 

S = T2^ Im ""'- fe(m) < 3 > 

where the electron mass has been taken to be zero. Putting d 4 p = it dM 2 p^dpx dy , 
where M, pt and y are the invariant mass, transverse momentum and rapidity of the 
lepton pair respectively, this is written as 

' IR " ' - M ^ osh ^ T W«. (4) 



dM 2 p T dp T dy 12vr 3 M 2 



in the Boltzmann approximation. 

The information about the specific processes occurring in the thermal medium which 
produces the dileptons resides in the current correlator . In the QGP, the domi- 
nant channel for dileptons production is quark-antiquark annihilation. The rate for this 
process is obtained from the lowest order diagram contributing to the current correlator 
W^ v . In the limit of vanishing quark masses, the rate for the process qq — ► e + e~ is 
obtained as, 

'IR a e -M T coshy/T /g\ 



dM 2 p T dprdy 9 Air 3 

Now let us consider dilepton production in the hadronic medium. The process 
7r + 7T - — > e + e~ is known to be the most dominant source of dilepton production from 
hadronic matter. In order to obtain the rate of dilepton emission due to pion annihi- 
lation from Eq. (||) the electromagnetic current correlator is expressed in terms of the 
effective propagator of the vector particle in the thermal medium using vector meson 
dominance (VMD) so as to obtain 



r 3e 2 m* 4 



W q 2 



Imll^ 



(M 2 - mf) 2 + (Imn«] 2 j 



(6) 



where Imll^ is the self energy of the p arising out of interaction with excitations in the 
medium and m* 2 = m 2 + (Rell^) 2 . For a p meson propagating with energy po and three 
momentum p the width is given by 



2T 

1 + 



x In 



l-exp[-f(p + ^( 5 )^ 2 - S )] 
1 _ exp[-f(p -W(s)JpJ^~s)] 



(7) 
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where s = M 2 = p\ — p 2 and W(s) = — 4m^/s. It can be easily shown that this 
expression reduces to the usual expression for the rho width in the limit p — > 0, 



■7r(P0) 



.9, 



p 7T 7T 

48vr 



Po^ 3 (Po) 



/be(-j)/bb(-^ 



(8) 

with / BB (aO = [e* - l]" 1 . 

We had shown earlier [29| that the width of the u meson can be large due to various 
reactions occurring in the thermal bath. The most dominant process, among others, 
is uit — ► 7T7T [28]. In view of this we also take into account the in-medium decay 
u> — > e + e~. Weldon had shown that the production of an off-shell vector meson 
(V) of four momentum p (where p 2 = M 2 ) in a thermal medium and its subsequent 
decay into a lepton pair leads to the dilepton emission rate 



dR 
~dfp 



(2tt) 3 

As before this reduces to 

dR 3 



V->l+l~ 



1 



Imlip 



vr (p2 - m* 2 ) 2 + (Imn«) 2 



dM 2 pt dpT dy 47r 3 
with the vacuum decay width 



1 V^l+l 



M 2 Y V {M) 



1 V^l+l 



(M 2 - m* 2 ) 2 + M 2 Tl{M) 
given by 

Ana 2 my*4 



-Mt coshy/T 



M 3 



(9) 



(10) 



(11) 



4 Space-Time Evolution 

The observed spectrum originating from an expanding hadronic matter is obtained by 
convoluting the static rates given above with the expansion dynamics. The usual space- 
time picture of the evolution is as follows. Matter initially formed as a hot hadronic 
gas expands and cools till freeze-out. However, in a phase transition scenario, matter 
produced as QGP expands and cools to the phase transition temperature T c . A mixed 
phase of coexisting quark and hadronic matter follows. After all the quark matter has 
converted to hadrons, the hot hadron matter expands and cools till freeze-out. In this 
work we use Bjorken-like |3lJ] hydrodynamical model for the isentropic expansion of the 
matter. The essential input at this stage is the equation of state (EOS) which provides 
the cooling law. For the QGP sector we use the bag model equation of state with two 
flavour degrees of freedom. The temperature in the QGP phase evolves according to 
Bjorken scaling law T 3 r = Tf ti. 

The hadronic phase is taken to consist of tt, p, u, r] and a\ mesons and nucleons. 
The energy density and pressure for such a system is given by, 



i=mesons 



(2tt) 2 



d 3 P E i f BE (E i ,T) + 



d 3 pE N f FD (E N ,T) (12) 



and 



■t— mesons 



(27T) 



,. / a Pjj^fBE{Ei,T) + 



(2k) 



, , d 3 p^-f F1) {t\.T) (!.",) 
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where t he sum i s over all the mesons under consideration and N stands for nucleons and 
Ei = \Jp 2 + mf. The entropy density is parametrized as, 

s H = €H+ T PH = 4a eff (T)T 3 = 4^ 9eff (m*,T)r 3 (14) 

where g e s is the effective statistical degeneracy. Thus, we can visualize the finite mass of 
the hadrons having an effective degeneracy g e R (m* , T) . (For a possible parametrization 
of entropy density consistent with lattice QCD we refer to p^| ). We will see later 
that the value of g e s is large at high temperature. The variation of temperature from 
its initial value Tj to final value Tf (freeze-out temperature) with proper time (r) is 
governed by the entropy conservation 

s(T)r = siTijn (15) 

Similarly, the evolution of the baryon density is governed by the conservation of baryonic 
current. In the present scenario the velocity of sound, c s becomes a function of T, 



_ 2 T dsn 



+ 3 



5eff dT 



(16) 



We will see later that c s differs substantially from its value corresponding to an ideal 
gas. We have considered a first order phase transition scenario where the the entropy 
density in the mixed phase is given by 

Smix = /q sq(T c ) + (1 - fo) s H {T c ) (17) 

where sq (s#)is the entropy density in the QGP (hadronic) phase at temperature T c and 
the function /q (1 — /q) is the volume fraction of QGP (hadronic) phase, obtained by 
solving the equation which governs the isentropic expansion of the matter. The initial 
temperature of the system is obtained by solving the following equation 

diV 7r _45C(3)_ D2 ,„ ^ 3 



dy 2vr 4 



n R 2 A 4a eS T? Ti (18) 



where dN n /dy is the total pion multiplicity (~ 1.5 x charge multiplicity), Ra is the 
radius of the system, r« is the initial thermalization time and a e g = (7r 2 /90) g c d (Tj). 
When the system is produced in the QGP phase, g e g in the above equation is replaced 
by gQGP which for two quark flavours is 37. If the quark and gluon masses are non-zero 
in the thermal bath then the effective degeneracy of the QGP, 5qqp, defined by Eq. (14) 



can be lower than 37, resulting in a higher value of Tj for a given multiplicity according 



to Eq. (18). Note that the change in the expansion dynamics as well as the value of 
the initial temperature due to medium effects relevant for a hot hadronic gas also enters 
through the effective statistical degeneracy. 

5 Results 

We will now compare our results with experimental spectra obtained by the CERES 
collaboration in Pb+Au collisions at the CERN SPS. In the central rapidity region the 
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entropy per baryon is found to be quite large ~ 40 — 50 [33, 34], consequently the 
effects of finite baryon density will be small. We have confirmed this by assuming an 
initial baryon density of n l B = 2.5n B for dN^/dr] = 150, 210, 270, and n B = A.5n B 
for dN c h/drj = 350, where n B is the baryon density of normal nuclear matter. The 
baryonic contribution is further suppressed due to the limited acceptance of the CERES 
detector |35], |3(| (see also ]37| ) . The dilepton spectra for various charge multiplicities are 
evaluated both with and without a phase transition. In all the cases we will assume 
an initial formation time Tj = 1 fm/c and freeze-out temperature, Tf = 120 MeV [38]. 



The initial temperatures depend on the scenario concerned. They are listed below in 
Table 1. The second column gives the initial temperatures in a phase transition scenario 
with QGP as the initial state. The next two columns are the values corresponding to 
a hot hadronic gas scenario with vacuum and QHD masses respectively. When mass 
modifications are incorporated using universal scaling the initial state depends on the 
particular value of T c that one considers. In view of the prevailing uncertainties in its 
value we have considered T c =170, 190 and 200 MeV @. For dN ch /dr] = 150 and 210, 
Tj < 170 MeV. Therefore, for these cases we consider the hadronic initial states only. 
However, for events with dN c f l /drj=270 and 350 the initial temperature can be larger 
than the critical temperature as seen below. 





QGP 


vacuum 
mass 


QHD 


Scaling 
T c =170 MeV 


Scaling 
T c =200 MeV 


dNch/dr] 


Tt (MeV) 


T t (MeV) 


Tt (MeV) 


Ti (MeV) 


Ti (MeV) 


150 


161 


209 


199 


168 


181 


210 


171 


219 


206 


171 


187 


270 


179 


227 


212 


179 


192 


350 


187 


235 


217 


187 


196 



18 self- 



Table 1 : Values of initial temperature obtained by solving Eqs. |14j and 
consistently for various values of dN c h/drj. 

The velocity of sound, which governs the space-time evolution of the system, is 
obtained by solving Eqs. (14) and (16) self consistently. This is shown in Fig. |l| where 
hadronic masses vary as a function of temperature according to universal scaling and 
QHD model. It is noted that the value of c 2 s in the hadronic phase is less than that of a 
massless ideal gas (c* rfea ' = l/\/3) indicating substantial interaction among the various 
constituents of the system. This is consistent with the results of Ref. [32|. 

In order to substantiate the point further, we now show in Fig. ^ how the effective 
degeneracy of the hadronic matter varies with temperature. The quantity g c s as obtained 
from Eq. ( [14] ) contains the effect of interactions among the various hadronic constituents 
which gives rise to the medium modified masses. We observe that when the masses vary 
according to Nambu scaling law, the value of g e g nearly approaches 30 which is more 
than thrice the value one obtains when vacuum masses are considered. As emphasized 
earlier, this variation has non-trivial effects on the cooling law of the system described 



by Eq. (|l5|) and is displayed in Fig. 
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Figure 1: Velocity of sound as a function of temperature, for different mass variation 
scenario for the hadrons in a deconfinement phase transition scenario. 
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Figure 2: The effective degeneracy as a function of temperature 
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Figure 3: The temperature as a function of time 
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Figure 4: Dilepton spectra for (A r c / l )=150. The dashed line (dash-dotted) indicates 
pair yield when the system is formed in the hadronic phase and hadronic masses 



e ' e 



vary according to universal scaling (QHD). Note that these also include the yield from 
hadronic decays at freezeout (dotted line). 



We have considered dileptons with transverse momentum above 200 MeV/c and 
an opening angle ee > 35 mrad. These are kinematical cuts relevant for the CERES 
detector and are incorporated as described in In all the figures, the quantity 

(iV c /j) indicates the average number of charged particles per unit rapidity interval in the 
pseudorapidity range 2.1 < r/ < 2.65. In all the cases background due to hadron decays 
are added to the thermal yields. 

In Figs. H and [B] we have shown the invariant mass spectra of dileptons corresponding 
to (N c h)=150 and 210 respectively. In these cases the initial temperature is not high 
enough to favour a phase transition to QGP. Hence we have considered a hot hadronic 
gas as the initial state. The medium effects are considered according to universal scaling 
(Eq. (||) with A = 1/2) and the QHD model. For dN^/dr] = 150 a good description 
of the data is realized with Tj = 168 MeV obtained by solving Eq. (18) for the mass 
variation given in Eq. (jl|) with A = 1/2 (dashed line). For QHD however, T { = 199 MeV 
(because of lower g e g) resulting in a larger dilepton yield. The theoretical yield (dash- 
dotted line) overestimate the data in the invariant mass region, M ~ 0.5 GeV because 
of smaller reduction in p-mass in QHD. For dN^/dr) = 210 the initial temperatures are 
Ti = 171 and 206 MeV for universal scaling and QHD model respectively. The dilepton 
yield is well reproduced for both the cases. Considering the experimental error bars it 
is not possible to state which one of the two is compatible with the data. 

In Fig. H we display the results for dN c h/drj = 270. In this case apart from the 
hadronic gas scenario we have also considered a deconfined initial state with Tj = 179 
MeV which evolves into a hadronic gas via a mixed phase. The observed enhancement 
of the dilepton yield around M ~ 0.3 — 0.6 GeV can be reproduced with the QGP initial 
state, once the variation of vector meson masses in the mixed and the hadronic phases 
are taken into account (solid line). The data is also reproduced by a hadronic initial 
state with Tj = 212 and 192 MeV for QHD (dash-dotted) and universal mass variation 
scenario (dashed) respectively. The p-peak in the dilepton spectra shifts towards lower 
M for universal scaling compared to QHD model, indicating a strong medium effect in 
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Figure 5: Same as Fig. | for (N ch )=210. 



the former case. In this case again it is impossible to differentiate between the QHD 
model and the universal scaling scenario. 

In Fig. [?] we depict the dilepton spectra for QGP and hadronic initial states for 
dN c h/dr] = 350. Results with hadronic initial state and universal scaling of hadronic 
masses with temperature seem to describe the data reasonably well. We have also noted 
that with the temperature dependent mass from QHD model can not reproduce the 
low M enhancement of the experimental yield. For T, = 187 MeV (see table 1) with 
a QGP initial state it is not possible to explain the measured dilepton spectra (result 
not shown in the figure). However, as we mentioned in section 4, the value of Tj in 
the 2-flavour QGP may be larger if gQGP is lower than 37 due to interactions among 
quarks and gluons (leading to non-zero masses of quarks and gluons) in the QGP phase. 
A good description of the data can be obtained by taking T, = 200 MeV with QGP 
initial state for T c = 190 MeV. All these results seem to indicate an initial temperature 
~ 200 MeV, a value which we obtained earlier by analyzing WA98 photon data. We 
also show the results due to large broadening of the p spectral function in the medium. 
The broadening of p is modelled by assuming the temperature dependent width as : 

r p (r) = r p (o)/(i-r 2 /r c 2 ). 

Finally, in Fig. || we demonstrate the effects of a large broadening of the vector 
meson spectral function on the photon spectra. Here we are interested in the relative 
shift of the photon spectra with vacuum mass and decay width of p in comparison to 
the spectra obtained with vacuum masses but large broadening of p. To estimate the 
photon yield from hadronic matter, we have considered the reactions [fiC| ], 7r p — > ir'y, 
it it —* p"f, ir ir —* 77 -7, ttt] — > 7r 7 and the decays p —* irirj and to — > ir'j. The 
invariant amplitudes for all these processes can be found in Refs. [25, 26]. We have 
also considered photon emission due to the reaction ir p — > a\ — > ttj [16]. The width 
of p has been taken as 1 GeV, independent of temperature purely for the purpose of 
demonstration. In a thermal medium, even if a huge enhancement of the width takes 
place as a result of interactions, near the freeze-out temperature the vector meson will 
regain its vacuum properties. So, by taking a large constant width throughout the 
evolution process we have overestimated its effect on the photon spectra. However, it is 
clear from Fig. ||, the effect is still small; the photon data with the current statistics can 
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Figure 6: Dilepton spectra for (N c / % )=270 for different initial states and mass variation 
scenarios. Solid line indicates dilepton yield when QGP is formed. (Here QGP scenario 
indicates sum of yields from QGP phase+QGP part of the mixed phase +hadronic part 
of the mixed phase +hadronic phase) The dashed line (dash-dotted) indicates e + e~ pair 
yield when the system is formed in the hadronic phase and vector meson masses vary 
according to universal scaling (QHD). 



not detect such effects. 

Throughout this work we have assumed thermal equilibrium of the system. How 
good is this assumption? We recall that the number of tt — ir collisions, N 7rn in the 



region M ~ m* can be estimated from the following equation [41|: 



N " = 2 x 3 x w„^-^ (19) 

where B p is the branching ratio for the decay p — > e + e~ . The factor 3 comes from the 
isospin combinations and counting collisions per particle gives a further factor of 2. We 
assume that thermal equilibrium may be realized in the system if N nn > 1 (a similar 
argument was put forward in ref. |^2[ regarding the equilibration of gluons). Considering 
the total multiplicity = 1.5 x dN c h/dr], we point out that thermal equilibrium may not be 
realized in the case of dN c h/dr] = 150 and 210 since the number of collisions per particle 
in the system, ~ N n7V / '(1.5dN c h/ 'drj) < 1. Interestingly, we note that for dN c h/dr] = 150 
the thermal source is not required to describe the data; the yield from hadron decays 
explain the data well. For dN c h/dr] = 210, however, the yield from hadronic decays 
underestimate the data. For higher multiplicities {dN c h/dr] = 270 and 350) the number 
of collisions per particle ~ 2 — 3. This is because for higher multiplicities (< > 
=270, 350) the pion density (~ T 3 ) increases by a factor of 1.6 and the number of it — it 
collisions is proportional to the square of the pion density. This indicates that thermal 
equilibrium may have been achieved, justifying the use of hydrodynamics to describe 
the space time evolution of the system. 

We have also checked that the condition, r s ^ att < r exp is satisfied throughout the 
evolution process, indicating that the use of hydrodynamics is reasonable for space-time 
description |)43| . Here Tg Cat t = (J2i( v ij) cr ij n j)~ 1 > i s the mean collision time, Vij is the 
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Figure 7: Dilepton spectra for (iV c /j)=350 for different initial states and mass variation 
scenarios. Solid line indicates dilepton yield when QGP is formed with T{ = 200 MeV 
and T c = 190 MeV. The dash-dotted line indicates e + e~ pair yield with p broadening 
and dashed line represents results when the system is formed in the hadronic phase and 
vector meson masses vary according to universal scaling. The dotted line shows dilepton 
yield from hadronic decays. 



relative velocity, cr^- is the scattering cross section, rij is the particle density of specie j 
and T cxp = TiiTi/T) 1 / " is the expansion time scale. 

It was proposed earlier [44] that the onset of deconfinement can be inferred through 
the dilepton yield if the quantity, dN ee /drj/(dN c i l /dr])' 2 (a quantity which is free from 
uncertainties in initial conditions) reaches a plateau with increasing dN c h/drj and the 
height of the plateau will be a measure of T c . We evaluated the above quantity for the 
four charge multiplicities mentioned above and did not observe any plateau. 



6 Summary and Discussions 

We have studied the dilepton yield measured by CERES experiment for various values of 
charge multiplicities in Pb + Au interactions. It is observed that for lower multiplicities 
{dN c h/dr] = 150 and 210), dileptons seem to originate from a hadronic source with initial 
temperature ~ 170 MeV. For the higher values of the charge multiplicity (270 and 350), 
however, the data can be described by both QGP and hadronic initial states. The value of 
the initial temperature realized in the high multiplicity events is similar to that obtained 
from the analysis of WA98 photon spectra. We find that both the enhancement in the 
photon and dilepton yields measured by WA98 and CERES collaborations respectively 
are due to a thermal source with initial temperature ~ 200 MeV. However, it is difficult 
to state which one of the two phases, QGP or hadronic gas is realized in the initial 
state of matter produced after the collision. We observe that the reduction of vector 
meson masses in the QHD model fails to explain the low mass dilepton enhancement for 
dN c h/dr] = 350. We also point out that the number of ir — tt collisions extracted from the 
dilepton spectra near m p are reasonably high in order to maintain thermal equilibrium in 
the system. We conclude by noting that in order to reproduce the data for large charge 
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Figure 8: Thermal photon spectra for Tj = 200 MeV and Tj = 1 fm/c. Solid (dashed) 
line indicates photon yield with large broadening of the p spectral function (vacuum 
spectral function). 

multiplicities either a large broadening of p or a substantial reduction in vector meson 
masses is required. Though it is difficult to draw a firm conclusion at this point, these 
phenomena are closely related to deconfinement or chiral symmetry restoring transition. 
Indeed, in Ref. |^5| it has been shown that a parametrization which mimics the dilepton 
emission rate from qq annihilation reproduce the CERES Pb+Au data well. At this 
critical and interesting juncture we look forward to the experiments at RHIC where 
larger and hotter systems are likely to be produced. Along with quantitative gains in 
all the signals there is a distinct possibility of electromagnetic radiation from the QGP 
phase which will pave the way for an unambiguous conclusion regarding the formation 
of this novel form of matter. 
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